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ABSTRACT 


A  three-dimensional  method  is  developed  to  improve  the 
computation  of  the  drift  force  and  moment  for  small-waterplane- 
area,  twin-hull  (SWATH)  and  surface  ships  in  oblique  waves 
with  zero  forward  speed.  Numerical  results  have  been 
computed  for  three  ships:  SWATH  6A,  Stretched  SSP,  and  MCM 
5371.  The  results  for  Stretched  SSP  show  good  agreement  with 
experiment.  For  MCM  5371,  the  results  of  two-  and  three- 
dimensional  methods  are  almost  identical  to  each  other  and 
these  results  show  good  agreement  with  experiment  when  the 
wavelength  ratio  is  not  too  small.  Even  though  there  are 
no  test  data  available  for  SWATH  6A,  the  application  of  three- 
dimensional  theory  is  likely  to  improve  the  results  of  drift 
force  and  moment  for  SWATH  ships. 

ADMINISTRATIVE  INFORMATION 

This  work  was  performed  under  the  General  Hydrodynamics  Research  Program  admin¬ 
istered  by  the  David  W.  Taylor  Naval  Ship  Research  and  Development  Center  (DTNSRDC) , 
Ship  Performance  Department  and  was  authorized  by  the  Naval  Sea  Systems  Command,  Hull 
Research  and  Technology  Office.  Funding  was  provided  under  Program  Element  61153N, 
Task  Area  SR0230101,  and  Work  Unit  1562-201. 

INTRODUCTION 

The  study  of  steady  and  slowly  varying  forces  and  moments  on  a  ship  in  oblique 
waves  is  important  to  maneuvering  and  seakeeping  analyses,  particularly  when  mission 
requirements  include  station  keeping  or  very  low  speed  maneuvering.  Recent  examples 
of  Navy  ship  designs  which  include  such  requirements  are  T-AGOS  (a  small -waterplane- 

area,  twin-hull  or  SWATH  ship)  and  a  monohull  mine  countermeasures  (MCM)  ship. 

1* 

Maruo  has  derived  the  equations  for  drift  forces  in  terms  of  scattered  wave 

amplitudes  for  the  two-dimensional  problem  and  in  terms  of  a  Kochin  function  for  the 

three-dimensional  problem.  However,  he  has  provided  no  numerical  computation  of  his 

2 

analytic  results.  Newman  has  derived  the  equations  for  drift  forces  and  yaw  moment 

for  the  three-dimensional  problem  using  a  slender  body  assumption  and  computed  forces 

3 

and  moment  for  a  Series  60,  C  *  0.6,  Lee  and  Kim  have  recently  derived  the  equation 

O 

for  the  roll  moment  in  addition  to  drift  forces  and  yaw  moment  and  computed  forces 
and  moments  for  a  SWATH  ship  and  a  monohull  ship  through  a  two-dimensional  motion 
program. 

*A  complete  listing  of  references  is  given  on  page  23. 


The  present  study  consists  of  the  application  of  three-dimensional  motion 
results  to  linear  theory.  The  force  and  moment  equations  are  derived  with  appli¬ 
cation  of  the  rate  of  change  of  linear  momentum  and  angular  momentum  previously 

2 

presented  by  Newman.  The  disturbance  velocity  potential  is  computed  with  the 

three-dimensional  motion  program  for  zero  forward  speed  developed  by  Hong  and 
4 

Paulling.  For  numerical  computation  of  drift  forces  and  yaw  moment,  three  models 
have  been  selected:  SWATH  6A,  Stretched  SSP,  and  MCM  5371.  The  results  for 
Stretched  SSP  show  good  agreement  with  experiment.  For  MCM  5371,  the  results  of  two- 
and  three-dimensional  methods  are  almost  identical  to  each  other,  and  these  results 
show  good  agreement  with  experiment  when  the  wavelength  ratio  is  not  too  small. 

Even  though  there  are  no  test  data  available  for  SWATH  6A,  the  application  of  three- 
dimensional  theory  is  likely  to  improve  the  results  of  drift  force  and  moment  for 
SWATH  ships. 


THREE-DIMENSIONAL  SOLUTION  OF  THE  VELOCITY  POTENTIAL 
The  coordinate  system,  oxyz,  is  defined  to  be  fixed  to  the  ship  with  no  forward 
speed.  The  oz-axis  is  directed  vertically  upward,  and  the  ox-axis  is  positive 
towards  the  bow.  The  oxy-plane  is  the  undisturbed  free  surface.  It  is  assumed  that 
the  fluid  is  incompressible  and  inviscid,  and  the  flow  is  irrotational.  The  total 
velocity  potential  of  the  fluid  in  incoming  waves  in  the  presence  of  the  ship  is 
represented  by 


$(x,y,z,t)  =  Re[4> (x,y,z)e~l£0t:]  -  Re[ (4>1+4>B) e-it0t ]  (1) 

where  <J>j  Is  the  potential  of  the  incoming  wave  and  is  given  by 


tbj.  «-  exp[Kz+iKx  cos  B-iKy  sin  3]  (2) 

and  <p0  is  the  disturbance  velocity  potential  due  to  the  presence  of  the  ship.  In 

Equation  (2),  A  is  the  amplitude  of  the  incoming  wave,  0)  is  its  frequency,  g  is  the 

gravitational  acceleration,  3  is  the  angle  of  incoming  wave  relative  to  the  x-axis 

2 

(3  “  0,  following  sea  and  180,  head  sea),  and  K  *  oo  /g  is  the  wave  number.  The 
potential  consists  of  the  following  velocity  potentials  as 
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Lv\: 


j-1 


(3) 


is: 


s>3  ^ 


El 

Xx 


where  ip  (j-1, 2, . . . ,6)  is  the  velocity  potential  arising  from  the  motion  of  ship  with 
unit  amplitude  in  each  of  six  degrees  of  freedom,  and  (j=l,2, . . . ,6)  is  the  ampli¬ 
tude  of  motion  in  each  of  six  degrees  of  freedom.  The  diffraction  potential  J «a 
represented  by  <p-j. 

The  potential  is  determined  as  the  solution  of  the  following  conditions: 

1.  Laplace  equation  in  the  fluid  domain 


3x 


a2  *2 
<p  d  $  3  p 

— **-  + - ■*-  +  - J-  s  n 
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for  j-1,2, . . .  ,7  (4) 
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2.  the  body  boundary  condition 
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for  j-1,2, . . . ,6  (5) 


for  j-7 


(6) 
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K 


the  linearized  free-surface  condition 


dtp 


-3. 


=  0 


on  z=0 


(7) 


The  right-hand  side  of  Equation  (5)  is  the  normal  velocity  component  at  the  ship's 
surface  and  the  unit  normal  vector  is  directed  into  the  fluid  domain.  In  addition 
to  Equations  (5)  through  (7),  i p^  must  satisfy  the  radiation  condition  of  outgoing 
waves  at  infinity  and  become  zero  as  z  becomes  infinitely  negative. 

The  solution  of  Equation  (4)  with  the  boundary  conditions,  Equations  (5) 

4 

through  (7),  has  been  presented  by  Hong  and  Paulling  in  terms  of  Green's  function 


as: 


3 


The  integral  in  Equation  (8)  is  evaluated  over  the  average  wetted  surface  of  the 
body,  Sg.  The  unknown  source  density  is  Qj ,  and  G  is  the  Green  function  for  the 
problem.  The  Green  function  for  a  pulsating  source  in  deep  water  can  be  written  as 


presented  by  Wehausen  and  Laitone  as: 


G(x,y,z;x',y',z')  =  -jp  +  -g-  +  2K  pv 

1  2 


a 

i 


y(z+z' ) 

M-K  Jo(,Jr) 


+  127TK  eK(z+Z,)  J  (Kr) 


(9) 


where  J  is  the  Bessel  function  of  the  first  kind  of  zero  order  and 
o 


Rx2  =  (x-x')2  +  (y-y')2  +  (z-z')2 


(10) 


R^2  =  (x-x')2  +  (y-y')2  +  (z+z')‘ 


(ID 


r2  =  (x-x')2  +  (y-y') 


^2 


(12) 


The  source  density,  Qj  can  be  found  by  applying  the  body  boundary  conditions.  From 


Equations  (5)  and  (6),  the  source  density  may  be  represented  as: 


-2tt  Qj  (x,y,z)  +  -|^  JJ  Qj  (x'  »y'  ,z')  G(x,y,z;x',y',z')  dS 


■i“nj> 


for  j-1,2, . . . ,6  (13) 


.  .  .*•  V*  -s  /*V-  N' 

*/  v*  v*  v‘  v*  O  v'  O  O  ~  « 


for  j  =  7 


(14) 


=  • 


The  numerical  solution  of  these  equations  is  given  in  Reference  4. 


EQUATIONS  FOR  DRIFT  FORCES  AND  MOMENT 

The  surface  which  enclosed  the  fluid  domain  is  divided  into  three  subregions: 
Sg  ship's  wetted  surface,  Sp  free  surface  and  S  control  surface  at  infinity  (see 
Figure  1).  The  rate  of  change  of  the  linear  and  angular  momentum  is  given  by 
(Reference  1) : 


dM 
_ x 

dt 


JJ 


[pn  +pV  (V  -U  )] 
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Cp(xny-ynx)+p(xVy-yVx) (Vn~Un) ]  dS 


VSF+S" 


(15) 


(16) 


(17) 


where  p  *  fluid  pressure 

p  •  density 

(nx>ny,nz)  =  components  of  the  unit  vector  n 

U  =  normal  velocity  on  S_,  and 

n  B  + 

(Vx,Vy,Vz)  *  components  of  the  velocity  V 

On  the  body  surface  and  free  surface,  V  *  U  .  The  contribution  from  the  pressure 

n  n 

term  on  S_  is  the  vector  of  forces  and  moment  acting  on  the  body: 
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dS 


(18) 
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Figure  1  -  Control  Surface 
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®  on  ®  an^  P  *  0  on  Sj,,  the  forces  and  moment  become  as  follows 
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».  * 

the  time  averages  of  Equations  (21),  (22),  and  (23)  are  the  mean  drift  forces  and 
moment.  Since  the  time  average  of  linear  momentum  of  the  periodic  motion  is  zero, 
the  drift  forces  and  moment  are  given  as 


F 
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(24) 
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JJ  IPVPVyV  dS 
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(25) 


[p(xny-ynx)+p(xVy-yVx)Vn]  dS  (26) 

The  pressure  in  these  equations  is  given  by  Bernoulli’s  equation 

p  -  -  p  Re(-i6j4»e~iCOt)  -  y  p  |  V |  2  -  pgz  (27) 

The  control  surface  consists  of  a  circular  cylinder  of  large  radius  about  the  z- 
axis  and  extending  from  the  free  surface  to  z  *  -  ®.  On  S^,  the  boundary  conditions 
are  represented  in  cylindrical  polar  coordinates  with  x  =  r  cos  9  and  y  *  r  sin  9. 
Through  a  substitution  of  the  pressure  terms  and  taking  the  time  average  of  the 
velocity  potential  the  drift  forces  and  moment  are  given  by 

2  f 27T 

Fx  *  -  g^'  I  |H(^)|2  cos  9  d9  +  cos  3  Re[H(ir-8)]  (28) 
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F 
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(29) 


r2lT 

\  -  -  ||  Im  H*(tt+6)  H'  (ir+0)  d0  +  ^  Im[H’(ir-e)]  (30) 

J0 

The  details  of  derivation  are  given  in  Reference  1  and  in  Appendix  A.  Here,  H(0)  is 
the  Kochin  function  and  defined  as 

exp(Kz+iKx  cos  6  +  iKy  sin  0)  dS  (31) 

The  asterisk  (in  Equation  (30))  denotes  the  complex  conjugate,  the  prime  denotes  the 
derivative  with  respect  to  0,  and  Re  and  Im  mean  the  real  and  imaginary  parts  of  the 
complex  value. 

If  we  substitute  Equation  (3)  for  <j>g,  the  Kochin  function  can  be  written 

6 

H(0)  -  H?(0)  Hj  (0)  (32) 

where 


H<9)  -  -  J  JGr  fe) 


Hj(0)  ■  -  J*J*  -  V’j  exp[Kz+iK(x  cos  0+y  sin  0)]  dS 


The  exciting  forces  can  be  expressed  in  terms  of  the  Kochin  function  (see 
Appendix  B)  as 


Fj  =  ^  [-Im  H  (-0)  +  i  Re  H  (-0)] 


Equation  (34)  is  useful  in  order  to  check  the  numerical  evaluation  of  Kochin 
functions  by  comparing  with  the  motion  results. 
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RESULTS  AND  DISCUSSIONS 


In  order  to  validate  the  numerical  results  of  the  derived  equations,  three  hull 
forms  have  been  selected:  SWATH  6A,  Stretched  SSP,  and  MCM  5371.  The  last  hull 
form  is  a  surface  ship.  The  principal  parameter  values  of  these  three  hull  forms 
are  given  in  Table  1. 


TABLE  1  -  VALUES  OF  PRINCIPAL  PARAMETERS* 


Parameter  (and  Unit) 

SWATH  6A 

Stretched  SSP 

MCM  5371 

Displacement  (long  ton) 

2579 

618 

1725 

Characteristic  Length  L  (m) 

54.3 

45.5 

73.2 

Length  of  Waterline  (m) 

52.5 

42.4 

— 

Length  of  Main  Hull  (m) 

73.2 

46.3 

— 

Beam  of  Each  Hull  at  Waterline  (m) 

2.2 

3.4 

13.4 

Draft  at  Midship  (m) 

8.1 

5.0 

3.6 

Maximum  Diameter  of  Main  Hull  (m) 

4.6 

3.4 

— 

Hull  Spacing  (m) 

22.9 

12.2 

— 

Longitudinal  Center  of  Gravity 

Aft  of  Main  Hull  Nose  (m) 

35.5 

21.9 

37.9 

Longitudinal  Center  of  Floatation 

Aft  of  Main  Hull  Nose  (m) 

35.2 

24.6 

Vertical  Center  of  Gravity  (m) 

10.4 

4.9 

Longitudinal  GM  (m) 

6.8 

18.7 

Transverse  GM  (m) 

3.4 

1.5 

Radius  of  Gyration  for  Pitch  (m) 

17.1 

12.4 

Radius  of  Gyration  for  Roll  (m) 

12.1 

5.6 

2 

Waterplane  Area  (m  ) 

193.9 

68.4 

739.0 

Length  of  Strut  (m) 

52.4 

31.5 

— 

Strut  Gap  (m) 

— 

10.8 

— 

Maximum  Strut  Thickness  (m) 

2.2 

1.2 

— 

♦Dimensions  are  full  scale. 

The  numerical  results  of  SWATH  6A  are  shown  in  Figures  2-5.  Figures  2  and  3 
show  the  results  of  head  sea;  Figures  4  and  5  show  the  results  of  beam  sea.  At 
head  sea,  either  motion  results  or  results  of  sway  drift  force  show  good  agreement 
between  two-  and  three-dimensional  results  except  for  small  X/L  values  where  the 
two-dimensional  theory  diverges  from  its  asymtotic  value.  In  beam  waves,  the 
results  of  sway  drift  force  and  yaw  drift  moment  show  some  discrepancies  between  two 
and  three-dimensional  approaches.  Without  the  experimental  results,  it  is  impossibl 
to  judge  which  one  is  correct.  These  discrepancies  occur  even  in  the  results  of 
restrained  motion  (if  this  is  the  case,  only  the  diffraction  potential  is  included 
in  the  computation).  This  indicates  that  in  the  solution  of  the  diffraction  poten¬ 
tial,  the  longitudinal  interaction  of  the  sections  plays  an  important  role.  The 
strut  length  of  SWATH  6A  is  52.4  m  and  the  length  of  the  main  hull  is  73.2  m.  Where 
the  strut  ends,  one  section  is  fully  submerged  and  the  other  is  floating  in  the  two- 
dimensional  model.  The  interaction  of  these  two  sections  is  not  included  in  the  two 
dimensional  approach;  however,  in  the  present  method,  this  interaction  is  fully 
taken  into  account.  The  full  effect  of  this  interaction  will  be  shown  in  a  follow¬ 
ing  example  for  the  Stretched  SSP.  The  two-  and  three-dimensional  results  for  heave 
and  pitch  motions  show  good  agreement  with  each  other  and  these  results  for  heave 
motion  are  in  good  agreement  with  test  data.  However,  the  analytical  results  for 
roll  motion  show  large  discrepancies  from  the  test  results.  These  discrepancies  are 
probably  due  to  the  effect  of  fins  since,  in  the  analytical  method,  the  effect  of 
fins  is  not  included  in  the  motion  results. 

Figures  6-7  show  the  computational  and  experimental  results  for  the  Stretched 
SSP  in  beam  waves.  The  computed  sway  forces  here  are  compared  with  the  experimental 
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results  presented  by  Numata.  There  are  some  minor  differences  between  the  three- 
dimensional  results  and  test  results;  however,  the  general  trend  is  the  same.  The 
two-dimensional  results  show  major  differences.  As  shown  in  the  motion  results, 
especially  for  roll  motion,  the  two-dimensional  results  show  large  discrepancies 
from  the  three-dimensional  results  and  the  test  data.  The  agreement  of  roll  motion 
between  three-dimensional  results  and  the  test  is  very  good;  however,  there  are  some 
discrepancies  in  heave  motion  when  X/L  is  larger  than  1.5.  These  discrepancies  do 
not  affect  the  results  for  drift  forces  since,  when  X/L  is  larger  than  1.5,  the  sway 
drift  forces  become  very  small  for  the  three-dimensional  case  and  test  data.  This 
model  has  a  large  gap  in  the  strut  in  the  vicinity  of  the  midship  section,  as  shown 
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Figure  3  -  Surge  Drift  Force  of  SWATH  6A  In  Head  Sea  at 
Zero  Forward  Speed 
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Figure  5  -  Sway  Drift  Force  and  Yaw  Drift  Moment  of  SWATH  6A 
in  Beam  Sea  at  Zero  Forward  Speed 
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Figure  7  -  Sway  Drift  Force  of  Stretched  SSP  in  Bern  Sea  at 

Zero  Forward  Speed 


in  Figure  8.  In  the  two-dimensional  modeling,  the  sections  which  have  struts  are 
modeled  as  floating;  whereas  those  sections  in  the  gap  region  are  fully  submerged. 
These  interactions  are  not  considered  in  the  two-dimensional  modeling.  The  three- 
dimensional  approach  takes  these  into  account  and  models  the  interaction  of  sections 
fully.  The  discrepancies  in  roll  motion  between  the  two-  and  three-dimensional 
approaches  are  probably  the  result  of  these  modeling  differences  and  the  interactions 
of  these  different  sections.  The  results  for  SWATH  6A  snd  Stretched  SSP  indicate 
that  the  sectional  interaction  significantly  affects  the  motion  and  drift  forv.e 
for  the  complicated  shapes  of  these  hull  forms. 

Figures  9  and  10  show  the  results  of  MCM  5371  in  beam  sea.  Even  though  there 
are  some  minor  differences  in  predicted  roll  motion  between  the  two-  and  three- 
dimensional  methods,  the  drift  forces  are  nearly  identical.  Contrary  to  the  SWATH 
ships,  this  hull  has  all  floating  sections.  When  the  sections  are  of  one  type,  in 
this  case  floating  (surface  intersecting),  the  interaction  of  the  sections  does  not 
seem  to  affect  the  motion  and  drift  force  results.  Compared  with  the  experiment 
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presented  by  O'Dea,  the  computed  results  of  drift  force  show  good  agreement  when 
X/L  is  larger  than  0.3.  As  X/L  values  become  smaller,  the  discrepancy  between  the 
present  method  and  experiment  increases.  The  heave  motion  results  show  good  agree¬ 
ment  with  the  test  results.  However,  the  measured  roll  motion  is  too  high  and  these 
are  not  plotted  in  Figure  9. 

SUMMARY  AND  CONCLUSIONS 

The  three-dimensional  method  is  applied  to  improve  the  computation  of  drift 
force  and  yaw  moment  for  two  SWATH  ships  and  a  surface  ship.  From  the  present 
study,  the  following  conclusions  may  be  drawn: 

1.  The  numerical  results  of  the  three-dimensional  method  for  the  Stretched  SSP 
show  good  agreement  with  experiment,  whereas  the  result  of  applying  the  two- 
dimensional  method  is  unsatisfactory  for  this  model.  Although  the  test  results  for 
SWATH  6A  are  not  available  for  validation,  it  appears  that  for  general  SWATH  ships 
the  application  of  three-dimensional  theory  is  more  realistic,  and  the  method  shows 
better  agreement  with  experimental  data  than  the  two-dimensional  method. 

2.  For  the  surface  ships,  the  results  of  two-  and  three-dimensional  methods 
are  almost  identical  and  these  results  show  good  agreement  with  test  data  when  X/L 
is  not  too  small. 


SWATH  6A  WET  DECK 


Figure  8  -  Side  Views  of  SWATH  6A  and  Stretched  SSP 


Figure  9  -  Heave  and  Roll  Motions  of  MCM  5371  in  Beam  Sea  at 

Zero  Forward  Speed 
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3.  At  the  present  time,  the  computational  time  for  the  three-dimensional 
program  is  fairly  high.  The  computational  time  can  be  reduced  in  the  future  by 
reorganizing  the  present  program  and  including  creation  of  a  master  table  of  the 
Green  functions. 

4.  The  present  study  is  restricted  to  the  computation  of  drift  force  and 
moment  for  zero  forward  speed.  The  investigation  of  nonzero  forward  speed  is  highly 
recommended . 
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APPENDIX  A 

DERIVATION  OF  EQUATIONS  FOR  DRIFT  FORCES  AND  MOMENT 


The  derivation  of  Equations  (28)-(30)  is  given  in  Reference  1.  Since  the 
potential  of  the  incoming  wave  is  expressed  differently  in  this  presentation  from 
that  given  in  Reference  1,  more  details  of  derivation  are  presented  in  this  appendix. 
In  the  following,  an  asterisk  denotes  the  complex  conjugate,  the  prime  denote  a 
derivative  with  respect  to  0,  and  Re  and  Im  denote  the  real  and  imaginary  parts  of  a 
complex  function,  respectively. 

In  the  polar  coordinate  system,  the  velocity  components  are 


V  «  V  cos  9  -  Va  sin  0 
X  K  U 


(A.  1) 


V  ■  V  sin  0  +  VQ  cos  0 
y  R  0 


(A. 2) 


Vn  "  ‘  VR 


(A. 3) 


Substituting  these  equations  into  Equations  (24),  (25),  and  (26),  we  obtain 


Fx  ■  -  JJ  tp  cos  0+pVR(VRcos  9-VgSin  0)]  R  d0  dz 


(A. 4) 


Fy  *  -  jj  [p  sin  0+pVR(VRsin  8+VqCos  9)]  R  d0  dz 
S’ 


(A. 5) 


s*  ■  '  JJ 


P  R  d9  dz 

R  t) 


(A.  6) 


where  the  superscript  bar  denotes  time  averaging  over  the  wave  period.  The  radial 
and  tangential  velocity  components  are 


(I  '■l“t) 


(A. 7) 
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Re 


dip  — iCJt^ 


(r  30  6 


(A. 8) 


The  far-field  potentials  are 

exp[Kz+iKR  cos (8+6)] 


(A.  9) 


<t> 


B 


H(tt+0)  exp 


^Kz+iKR+i 


(A. 10) 


Equation  (A.  9)  Is  the  expression  of  Equation  (2)  in  polar  coordinates  and  is 
given  in  Reference  1.  Since,  in  the  far-field,  both  (j>T  and  4>_  are  exponential 
functions  in  z,  it  is  convenient  to  integrate  Equations  (A.4)-(A.6)  with  respect  to 
z.  The  upper  limit  of  this  integration  is  the  free  surface  height 


z  -  C  - 


-  -  * 
g  t 


Re  <)>(x,y,o)e“ia)tj 


z*o 


(A. 11) 


and  it  follows  that 


f  p  dz  -  p  C  Re[ico^(x,y,o)e-ia>t]  -  |  V |  2  -  y  pg£2  (A. 12) 

T-00 

The  first  and  the  third  terms  of  Equation  (A. 12)  become 


|pC  Re t iox}» (x,y,o)e”ia)t]-  l j 
-  (£  [Reac^e-1"')!2-  f  [Re  [f  *  e'1”')]  *  j 

‘  [Re(iaj^e*i(ut)  ]  *  ^  I  ^_0 
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The  second  term  becomes 


-«  w2  -  -  fe(|M|2+^  |H|2+k2w2)*-o 


(A. 14) 


The  Integration  of  and  V^Vg  in  Equations  (A.4)-(A.6)  are  given  by 


J  vr2  dz  ■  k  <v 


2,  .  1_  lit 

R  'z«o  4K  3R 


(A. 15) 


±  _ 

f  v  V  dl.L(vT)  -  JL  (ii  Ml  +  i*l 
J  R  8  2K  ^  R  0'z«o  8KR  V3R  39  3R  39/ 


(A. 16) 


Substituting  Equations  (A.13)-(A.16)  into  Equations  (A.4)-(A.6),  we  obtain  the 
forces  and  moment  in  terms  of  the  velocity  potential 


V  m  £ _ 

x  8K 


J  (^2  H  IF  -  H  IF  -r2<m>*)  R  COS  9  d6 


+  £- 
8K 


m  2fr 

(  (ii  Q)  8in  6  de 

J  \  3R  30  +  3R  30  /  810 


(A. 17) 


F  -  fn- 
y  8K 


2ir 

I  (? 


3^  1*  _  3£  3£*  2  \ 

3<|)  30  3R  3R  *  99  / 


R  sin  0  d0 


_  £_ 

8K 


/•27r 


(A. 18) 
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(A. 19) 


M  =  -  £— 
z  8K 


In  the  above  equations,  the  potential.  Its  derivatives,  and  their  conjugates  are 
computed  at  z  ■  0.  Employing  the  definition  <j>  =  <j>T  +  <J>„,  Equations  (A. 9)  and 
(A. 10)  are  substituted  Into  the  force  and  moment  equations.  The  first  term  of 
Equation  (A.  17)  becomes  zero  as  R  -*•  00 .  The  remainder  of  the  terms  can  be  expressed 
In  integral  form  as 


(•'"iim*  „„„  F  (**!  9**I  3**»  +  ffi  .  cos  e  d9 

J  9R  9R  R  COS  9  d0  "  J  \9R  9R  9R  9R  9R  9R  9R  9R  /  R  COS  9  d9 


*—  f 

2.  J 


H(ir+0)  H*(ir+6)  cos  0  d0 


/KR\1/2  f2lT 

-  2(t)  AK  (  Im  I  cos (0+0)  cos  0  H(ir+0)d0 


x  exp  £-1  KR  cos  (8+0)  +  1  KR+-1 
2tt  2tt 

I  K2  <p  <p*  R  cos  0  d0  ■  K2  I  (<J>j^j.*+<{ig<))g*+<)>j4,B*+4I*<l,g)  R  cos  0  d0 


(A. 20) 


—  f 

2tt  J 


H(tt+0)  H*(tt+0)  cos  0  d0 


1/2  r2Tr 

-  2oj  AK  Tm  I  H(ir+0)  cos  0  d0 


x  exp  [-1  KR  cos(8+0)  +-1  KR+i  f  1 


(A. 21) 


»2tt 


far  •*•«-/*'( 


3<J>T  34>*t  34>t  3<t>*  34>*t  9<J>*  3<J>*«  \ 

8R  55”  +  3R  §0  +  9R  W  +  3R  §5”  J8in  6  d6 


1/2  ^  ^ 

-  id)  AK  (H)  I  H(lT+0)  sin(6+e)  sin  9  exp  [-iKR  cos(B+0)+iKR+i  J  J  d0 


(A. 22) 


f2lT  3d>*  3d>  f2V^*T  3**T  3*R  3**r  3*T  3**B  *K\ 

J  H-  H  a  «*e  -  J  {- nr*r  +  ^rjr*-w-»r*-w-»r ) 


sin  6  d0 


-  i  (o  AK 


1/2 

(H)  J  H* (ir40 )  sin(8+6)  sin  0  exp  j^LKR  cos(B+6)-iKR-i  d0 


(A. 23) 


With  substitution  of  Equations  (A.20)-(A.23)  into  Equations  (A. 17)  and  (A. 18),  the 
forces  are 


F  - 

x 


.  £id 

8tt 


f2ir  1/2  r2lT 

I  |H(ir+0)  | 2  cos  0  d0  I®  H(it+0)  (cos  0+cos  8)  d0 


x  exp  £-IKR  cos(B+6)+iKR+i  jJ 


(A. 24) 


8tt 


/*2^  1/2  r 

j  jH(7r-h0) ) 2  sin  0  d0  +  (||)  Im  J 


2tt 


H(ir+0)  (sin  9-sin  B)  d0 


x  exp 


[' 


iKR  cos 


(B+0)+iKR+i 


(A. 25) 


The  integral  terms  of  Equation  (A. 19)  can  be  written  as 


f 2lT  M  Mi  r  ae  -  f27!!l^  +  ^^  +  ^^  +  ^^Rd9 

J  8R  96  R  d9  ■  J  ^gg-  9e  +  3r  ae  9R  36  +  8r  80  / 

0  o 

1/2  f2lT 

-  cuA  j  H* *  (tt+0)  cos(3+0)  exp  [lKR  cos(0+0)-iKR  -i  JJ  d0 

0 

1/2  f2lT  1 

-  ojA  (-  j  +1KR^  J  H(ir+0)  sin(0+0)  exp  £-iKR  cos(8+0)  +1KR+1  £jd0 


ijd  f 

2ir  J 


H(ir+0)  H*  ’  (tt+0  )  d0 


(A. 26) 


f2’  >t;  m  ,  „„  f 2”  /**i  8*i  +  **i !!» + !!!» !!i  + !!!» !!»)  *  d0 

I  3R  36  R  d6  *  J  \3R~  36"  +  TT-  36  +  3R  36  3R  39  / 

0 

1/2  r2lT  i 

.  ojA  (||)  J  H'  (tt+6)  cos(B46)  exp  [-1KR  cos(8+0)+iKR+i  Jjd0 

Jo 

1/2  2ir  1 
+  ojA  (H)  ^  +1Kr)  f  H*(ir+0)  sin(8+0)  exp  £  iKR  cos(8+0)-iKR-i  j  J  d0 

J0 

-<r 


h* (it+6 )  H'Cn+e)  de 


(A. 27) 


Substitution  of  these  two  equations  into  Equation  (A. 19)  yields  the  moment  as 


.2it 


f*”  1/2  r 

-  H  to  J  H'Or+6)  d0  Re  J 


M 


2tt 


H'(7r+6)  d0 


x  [l+-cos(8+9)  ]  exp  £-iKR  cos(8+9)+iKR+-i 


(A.  2! 


The  integrals  in  Equations  (A. 24),  (A. 25),  and  (A. 28)  involving  R  can  be  evaluated 
by  the  method  of  stationary  phase.  When  R  is  large,  the  following  integral  is 
integrated  as  (see  Reference  6) 


j  HZ, R) 


eiR*  (C)  d£  «  £  ^(ct 


r,R>  [ - — - 1  exp[iR.p(a  )+i  rl 

[Rk"(ar)  |  J  L  r“  4J 


♦27*<vr>  ^ r  6 

S  8  7% 


’ _ 6 

r|v>”  * 


‘“s’ 


1/3 


exp[iR«p(a  )] 
s 


/  +  sign  when  ^"( a^)  >  0 


<p'  (ar)  -  0,  v>"(“r)  *  0 


(  -  sign  when  ^"(ar)  <  0 


<fi'  (a  )  -  0,  ^"(a  )  -  0,  <P'"(a  )  *  0 

8  S  S 


Therefore,  the  exponential  integral  becomes 


I 


2ir 


f (0)  exp[-iKR  cos (8+9)]  d9 


[f(-8)e-1KR+1(,,M)+f(1r-6)e-KR-1(’,/4>]  CA.J. 


(A. 2' 


JJJ  IfV? WU  ^  -*1”' V 


W.W  5  '-  '■'■  '■•  «•  -•-  .1,  *  ■" 


0. 


g. 


i 


and  the  integrals  in  the  equations  of  the  forces  and  moment  are  evaluated  as 
follows 

f2* 

I  H(ir+0)  (cos  0+cos  0)  exp  £-iKR  cos(0+0)+iKR+i  ^-J  d0 
0 

1/2 

=  (-gj  2i  cos  0  H(n-0) 

.2-n 

I  H(it+0)  (sin  0-sin  0)  exp  £-iKR  cos  (0+0)+iKR+i  J-J  d0 
0 

/  ^  \  1/2 

“  ■  (kI/  2i  sin  s  H(tH3) 

f2n 

J  H'(tt+0)  [l+cos(0+0)  ]  exp  |^-iKR  cos(0+0)+iKR+i  d0 
0 

-  (I) m  21  »•  <*-« 


(A. 31, 


(A. 32) 


(A. 33) 


With  the  substitution  of  Equations  (A. 31)-(A. 33)  into  Equations  (A. 24),  (A. 25),  and 
(A. 28),  respectively,  the  equations  of  forces  and  moment  are  given  in  final  form  as 


_  £*t 

8ir 


I 


2-n 


|H(w+0)|: 


cos  0  d0  +  cos  0  Re[H(ir-0)] 


(A.  34) 


8tt 


I 


2-n 


|H(TT+0)  T  sin  0  d0  - 


sin  0  RefH(7r-0)]  (A. 35) 


M 


L. 

8tt 


I 


2-n 


H* (tt+0)  H’(tt+0)  d0  +  Im[H*  (tt— 0 )  ] 


(A. 36) 


APPENDIX  B 

DERIVATION  OF  EXCITING  FORCES  IN  TERMS  OF  KOCHIN  FUNCTION 


The  exciting  forces  are  computed  by  integrating  the  pressure  due  to  incoming 
wave  diffraction  potentials  over  the  wetted  surface  as 


Fj  -  -  jj  p  n^  dS  =  jj  (-ipco)  (^-K^)  n^  dS 


(B.l) 


From  the  body  boundary  condition  (Equation  (5)),  Fj  becomes 


Fj =  p  JJ  v  ( W  ds 


(B.2) 


Applying  Green's  identity,  we  obtain 


rj ' D  If  (vi^j  sir) 


(B.3) 


From  the  body  boundary  condition  for  the  diffraction  potential  (Equation  (6)),  F_.  is 
given  by 


Fj  •  pIf  (v*i  Is)  *i dS 


(B.4) 


We  substitute  Equation  (2)  for  and  finally  have  the  exciting  forces  as 


_  iPeA 


II  (vi  Is) 


e  exp(iKx  cos  g-iKy  sin  g)  dS 


^  Hj  (-B) 


(B.  5) 
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